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ABSTRACT 

Transit light curves for stellar continua have only one minimum and a "U" shape. By contrast, 
transit curves for optically thin chromospheric emission lines can have a "W" shape because of stellar 
limb-brightening. We calculate light curves for an optically thin shell of emission and fit these models 
to time-resolved observations of Siiv absorption by the planet HD209458b. We find that the best fit 
Si iv absorption model has R p ,siiv/R* = 0.34^o'?2? similar to the Roche lobe of the planet. While the 
large radius is only at the limit of statistical significance, we develop formulae applicable to transits 
of all optically thin chromospheric emission lines. 

Subject headings: planets and satellites: atmospheres — stars: chromospheres — ultraviolet: planetary 
systems 
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1. INTRODUCTION 

Since the first observation of a transitin g exoplanet 
(|Charbonneau et al. 2QQQ: lHenry et al.ll2000ft . knowledge 
of exoplanetary radii, composition and atmospheres has 
grown explosively. In-transit and out-of-transit spec- 
troscopy and photometry have revealed water absorp- 
tion in HD 189733b (Bea ulieu etHI 120 08). atmospheric 
emission in TrES-1 (iCharbonneau et al.ll2005ft. a surpris- 
ing n umber of anomalously large planets (Bar affe et al.l 
120101). and constraints on exoplanetary composition 
(Rog ers fc Seagerl 12010ft . Accurate transit timing can 
also indicate the presence of additional bodies in the sys- 
tem thxouj^^ or- 
bit ([Agol et al.ll2005l : iHolman fc Murravll2005ft . 

When transits are observed in stellar Lyman-a, Cn, 
Sim, Mgn and Oi emission, they show much deeper 
minima than for visible wavelengths, revealing escaping 
atmosphere s extending far beyond th e geometric radio 
of planets (iVidal-Madiar et al.l l200i iBen-Jaffell 1 2007 
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Murrav-Clav et alT 20091 ; iLecavelier des Etangs et aD 



2010; Fossati et al. 2010). These light curves also con- 
strain atmospheric c onditions and mass escape from the 
planet's Roche lobe dKnutson et aLll2007t iGarcfa Muhozl 
[20071 : iLinskv et al.ll2010ft . 

Transit searches have been largely in the optical and 
near-infrared continuum, where the star is optically thick 
and limb darkened due to the temperature profile at the 
r « 1 surface of the star. For limb darkened wavelengths, 
the flux from the transiting system is at a minimum 
when the planet crosses the sub-earth longitude of the 
star (phase = 0.0) because the stellar disk is brightest at 
its center. 
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By contrast, emission lines from stellar chromospheres 
and transition regions can be limb brightened. For op- 
tically thick emission, limb brightening occurs when the 
source function increases with stellar altitude. For op- 
tically thin emission, strong limb brightening occurs be- 
cause the chromospheric and transitio n region gas has its 
largest column density at the edges. Ass ef et all ([2009) 
showed that transit observations of chromospheric emis- 
sion lines decrease sharply to a minimum at the first limb, 
increases to a local maximum mid-transit and then re- 
verses the process as the planet exits the stellar disk. As 
a consequence , such a transit curve will be "W" -shaped. 

Ass ef et all (j2009ft point out that limb brightening 
could be useful for detecting exoplanetary transits of gi- 
ant stars. Light curves of limb brightened wavelengths 
have deeper minima than for both limb darkened and 
uniform disk emission. The star emits over a smaller ef- 
fective area-a ring instead of a disk-so the planet covers 
a larger amount of the stellar flux. This is important for 
transits of giant stars where broadband transit depths 
can be below 0.01% for Jupiter-sized planets. Also, the 
planet covers its host's limb for a small fraction of the 
transit, allowing for feasible detection of giant star tran- 
sits with ground based telescopes, which suffer from sys- 
tematic photometric errors over timescales longer than 
one night. The exoplanets 4UMa b, HD 122430b, HD 
13189b, and HIP75458 b all have transit probabilities 
greater than 10% and m ay be useful targets for future 
st udies (lAssef et al .112009ft . 

lAssef et al.l (j2009ft approximate the limb brightened 
star as a central disk of emission surrounded by a circu- 
larly symmetric ring with ^30 times the intensity. With 
this ring approximation, the maximum depth of the tran- 
sit is proportional to the ratio of the planet radius to stel- 
lar radius, R p /R* if the emission from the central disk 
is negligible, instead of (R p /R*) 2 , as expected for a uni- 
form disk. This is because the planet covers ~2R P out of 
a circle of emission whose total circumference is 2ttR*. 

In this paper, we present a transit light curve calcula- 
tion for an optically-thin and geometrically-thin shell of 
emission. The calculated maximum transit depth scales 
as (R p /R*) 3 / 2 , instead of R p /R*. In $2] we calculate 
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absorption profile is larger than the star. In this case, the 
same formulae apply-see Equation [2j A small figure for 
each regime is included in the equations, where an open 
circle represents the star and a filled circle represents the 
planet. 



S(x) = Q(p-x) + 



a 



(i) 



x +p 
x — p 



II(n,m) — AxpaiE(m) — a2K{m) 



Fig. 1. — Model Limb Brightening, (a) A model of spherically 
symmetric optically thin emission varies continuously from limb to 
center, (b) An approximate model where most of the emission is 
from the stellar limb surrounding a uniformly emitting circle. We 
employ the model shown in (a) for Si IV emission, whose transit 
light curve is given by equation [T] 



the expected limb brightened light curve for an optically 
thin emission line. We consider optically thin emission 
because it shows at least 8 times the limb bright ening of 
optically thick emission ([Kastner fc Bh atia 1992). We fit 
this model light curve to Siiv emission from HD209458 
in £[3] and discuss the implications for HD209458b's ther- 
mosphere in £14.11 

2. A LIMB BRIGHTENED CURVE UNDER THE 
THIN-SHELL APPROXIMATION 

We make the approximation that the thickness of the 
chromosphere, ft, is much smaller than the size of the 
planet (R p ) and star (i?*). Under this geometrically 
thin approximation, the total flux from the star is pro- 
portional to the surface area of the hemisphere facing 
the Earth times its thickness ft, because the total flux 
for an optically thin emission line is proportional to the 
total number of emitting ions. Neglecting any photo- 
spheric contribution, the amount of emission that the 
planet blocks is then simply the amount of the stellar 
surface that the planet covers times ft. In this geomet- 
rical limit, therefore, the thickness of the chromosphere, 
ft, cancels out. 

We compute the light curve for zero thickness (ft = 0) 
by finding the area of the the planet's shadow and divid- 
ing this by the surface area of a hemisphere with radius 
R*. The volume of the intersection of a cylin der with 
a sphere is given by lLamarche fc Leroyl (jl990f ) and we 
find the surface area of the intersection by taking a par- 
tial derivative with respect to the radius of the sphere. 
The result can be expressed analytically in terms of el- 
liptic integrals. Let x be the distance, in units of stellar 
radii, R* : from the center of the planet to the center 
of the star projected onto a plane perpendicular to the 
observer. Let p be the planet/star radius ratio, R p /R*. 
To calculate the light curve for a planet that does not 
pass through the center of the star, one can simply write 
x(t) as y/d 2 + b 2 , where d is the distance to the closest 
approach point in stellar radii R* and b is the impact 
parameter in stellar radii R*. 

The transit depth, S(x) is a piecewise function with 
three different regimes: (1) when the planet is fully 
contained in the stellar disk (2) when the planet is at 
egress/ingress and (3) when the planet is beyond the 
stellar disk. We also include the case that the planet's 
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where K(m), E(m), and H(n,m) are the complete 
Legendre elliptic integrals of the first, second, and 
third kinds, and Q(x) is the Heaviside step func- 
tion. For the elliptic integrals, we use the conven- 
tions of I Abramowitz fc Stegunl (jl972h where II(n,m) = 
U(n;K(m)\m) for the third elliptic integral. Q 

These formulae are difficult to evaluate numerically at 
x = 0, x = p and x = l±p due to the formal divergence 
of different terms in equations [T] and [2j the divergences 
cancel out analytically, but routines that evaluate the el- 
liptic integrals diverge. However, these locations are a set 
of measure zero, and thus are tractable when modeling 
data. 

2.1. Analytical Transit Depth Estimation 

The transit is deepest slightly before second contact, so 
we can estimate the maximum transit depth as follows 
from Figure [2] by comparing the total emitting area of 
the star to the total stellar surface area blocked by the 
planet. The blocked surface is the same as the shadow 
produced by a sphere of radius R p onto a hemisphere of 
radius R*. When the planet occults the edge of the star 
(second contact), as shown from an edge-on viewpoint in 
this figure, then the length of the arc of the long axis of 
the shadow is R*6. The diameter of the planet is 2R P « 
— cos#) & ^R*6 2 , where the latter approximation 
is valid for <C 1. We can approximate the shadow 
as an ellipse with a semi-minor axis of R p and a semi- 
major axis of 7}R*0, so the area of the shadow is A t = 
7T^R P R*R P . Thus, the maximum depth of transit is 

7 An IDL procedure for the transit depth 6(x) is located at 
http : //www . astro . Washington . edu/ agol/. 
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Fig. 2. — Edge-on view of area of stellar emission and the amount 
blocked by a planet. This blocked surface area is the same as the 
area of a shadow cast by a sphere onto a hemisphere. 



given by 



R p 



3/2 



(3) 



which is accurate to within 5% for R p /R*< 0.23. Note 
that this is a differe nt scaling fo r maxi mum transit depth 
than that given in lAssef et al.l ([2009) who assume that 
most of the stellar emission is from a thin ring. For 
Rp/R* = 1 — \/3/2, half of the stellar emission is within 
R p of the stellar radius i?*, so for R p /R* > 0.13, the 
ring approximation is valid, but the scaling of transit 
depth assumes negligible limb curvature at the scale of 
the planet (R p /R* 

The remarkable consequence of Equation [3] is that the 
depth of a chromospheric transit does not decline as 
much with the radius of the planet as a transit of a 
uniform disk. A chromospheric transit has a maximum 

/ \ 1/2 

depth that is « ^ (r) times deeper than the maxi- 
mum transit depth of a uniform disk; thus smaller plan- 
ets have an advantage to be obse rved at chromosp heric 
wavelengths, as emphasized by lAssef et all (|2009f) . It 
should also be noted that at mid-transit, the Double-U 
curve has a smaller transit depth than for a uniform disk, 
because the planet covers only ttR^ out of a hemisphere 
of area 2itR^. 

Figure 3 (solid curve) shows the transit light curve, 1- 
5(x), for a planet that has R p /R*=0.08. The estimate 
given by equation [3] and the more detailed equation [T] 
agree well, predicting maximum chromospheric depths 
to be 1.8 times deeper than for uniform disk brightness. 
(For a uniform disk emission 5 max = irR 2 p /{irRl)). We also 
include a light curve for the thin circle of emission shown 
in Figure [I] b for comparison. 

3. Si iv ABSORPTION BY HD 209458B 

iVidal-Madiar et al.l (|2003l 12004 ) observed the exo- 
planet host HD 209458 with the Hubble Space Telescope 
Imaging Spectrograph (STIS) instrument and found 
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Fig. 3. — Transit light curve for R p /R*= 0.08, using a family of 
3 different models. The solid line is for a transit of an optically 
thin shell using equation [T] wh ere the emission loo ks like Figure 
a. The "Uniform Emission" (Mandel & Agol 2002, dashed green 
line) is one for which the emission is assumed to be constant across 
the stellar disk. The "Thin Ring Emission" model (dot-dash red 
line) is a model where the emission is assumed to be mostly from 
a thin ring as pictured in Figure ^b. Note that the "Thin Ring" 
and "Optically Thin Emission" models' minima are deeper than 
the uniform emission model and that these minima occur near the 
stellar limbs. 




Fig. 4. — (a) SilV emission from the Solar limb as observed by 
SUMER (S olar Ultraviolet Measurement of Emitted Radiation) 
( Wiik et al. 1997) is strongly limb-brightened, indicating that a 
transit of this emission line should have a Double-U light curve, 
(b) Hell 304 A image of the solar limb, taken b y the Extreme- 
ultraviolet Imaging Telescope (Feldman et al. 2000) is, by contrast, 
not limb brightened because it is optically thick in the chromo- 
sphere and thus the additional column density at the limb does 
not contribute any more flux than the central disk. 



an extended hydrogen, oxygen and carbon atmosphere 
around the planet by fitting the lig ht curves to the Hi, 
O i an d C n lines during a transit. IVidal-Madjar et"aTl 
(2004) found no Siiv absorption when fitting the light 
curve to a model of a spherical planet occulting a uni- 
form stellar disk. 

We fit their Siiv 1394 A data (jVidal-Madiar et al.l 
(2004), Figure 3) with a Double-U model given by equa- 
tion Q] This model is appropriate for the Si iv emission, 
as evident in the Solar image i n Figure [4] (a) where strong 
limb-brightening is apparent. iFrance et al l (|2010l ) found 
that the Siiv 1394 A/ Siiv 1403 A ratio is 2:1 within the 
errors, ind icating that the Si iv 1394 A emission is opt i- 
cally thin (jBloomfield et al.ll2002 IChristian et al.ll2006D . 
We investigated the transit of the Siiv line because it 
the strongest optically thin line in HD209458b's STIS 
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spectrum. While other limb brightened emission lines do 
exist, we focus on the optically thin ones because they 
should be the most limb brightened. 

The model has two free parameters: the planet/star 
radius ratio R p /R* and an overall constant that sets 
the off-transit flux. The second parameter is necessary 
since the off-transit Si iv 1394 A STIS flux is poorly con- 
strained. The impact parameter is fixed with a value of 
b = 0.50i?*, using an inclination of 86.7 ° and a semi- 
major axis a = 8.76R* (jTorres et al.ll2QQ8[ ). 

Figure 5 show s the light curve for the 
IVidal-Madiar et al.l (|2004f ) data and the best fit 
model. In addition to the best fit model, we show 
two more limb brightened models where R p /R* is a 
fixed parameter for comparison. These have R p /R* 
equal to and 0.12, representin g no Siiv absorptio n 
and the geometric planet radius (Kn utson et al.l 12007). 
respectively. 

Using a Leyenberg-Marquardt fitting algorithm 
(jMarkwardtl l2009f ) and the points at which the 
X 2 — Xmin = 1 to calculate uncertainties, we find that 
Rp/R* = 0.34+Q12, close to the size of the Roche Lobe. 
The Double-U model has a total x 2 of 10-4 with 12 
data points and two parameters, which is 3.1 less than 
fitting the data to a constant flux line (representing 
no detection) with one parameter. Since these models 
are nested, we can use the maximum likelihood ratio 
test. The P- value for the difference in x 2 is 0.05, so the 
data favor the Double-U model with 95% confidence, 
assuming normally distributed data. The same test 
favors the R p /R* =0.34 model over the R p /R* = 0.12 
model with 91% confidence. 

4. CONCLUSION 

As indicated by the Solar image of Si iv in Figure HJ 
a limb brightened model should be used for Siiv emis- 
sion. The Double-U model for the Siiv transit fits the 
data better than a non-detection, even when accounting 
for the additional parameters in the model. The best 
fit absorption radius R p /R*=0.341q'i2? if the absorption 
is entirely optically thick. This radius favors planetary 
atmosphere models with mass flow beyond the planet's 
Roche lobe. 

The best time to observe transits in optically thin stel- 
lar emission lines is when a planet crosses its host's limb 
and not at a phase of 0, since the Double-U curve is deep- 
est at the stellar limb. Observations at the limb have a 
depth of ~0.5 (R p /R*) 3//2 whereas at the center of the 
star they have a transit depth of ~0.5 (R p /R*) 2 , only 
half the depth of a uniform brightness transit. 

4.1. Discussion 

The Siiv transit dept h is comparable to t h e O i tran- 
sit depth calculated by Vida l-Madjar et all ( 2004f) and 



larger than t he H i transit dep th calculated bv lBen-J affel 
(|2007D . As iKoskinen et all (|2010h point out, a hard 
sphere, which we assume in equation [TJ is a poor ap- 
proximation to planetary atmospheric absorption. The 
Siiv absorption profile, instead of having a sharp transi- 
tion from opaque to transparent, should change smoothly 
from optically thick to optically thin absorption in a more 
accurate model. In order to explain the large transit 
depth, this smooth model would have to have a radius 
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Fig. 5. — The light curve of the limb-brightened transition re- 
gion line Siiv during the transit of HD209458b, as taken by 
IVidal-Madiar et al.l (2004) for a wavelength range of [1391 A, 1397 
A] . Fluxes are normalized so that the weighted average of the points 
when the planet is not occulting the star is 1.0. The solid line is 
a best fit model, using equation [T] where the radius is a free pa- 
rameter. For reference, two different best-fit models are shown 
where the planet size is a fixed parameter: the dashed curve is a 
horizontal line, corresponding to R p /R* = or a non-detection of 
Siiv absorption. The dash triple-dotted line is a model with the 
geometric radius, R p /R* = 0.12. The best-fit radius, R p /R*=0.34 
(normalized x 2 = 1.07), indicates Siiv absorption from a Roche 
lobe-sized cloud. 



extending beyond the planet's Roche lobe to explain the 
obser ved transit depth. RR OC he/R* va ries from 0.35 to 
0.48 (jBen-Jaffel fc Sona Hosseffll2010h . 

The best-fit transit depth supports models with high 
concentrations of metallic ions in the atmosphere because 
the r adius of Si iv absorptio n is as large as for H i absorp- 
tion. IKoskinen et al.l (j201Q[ ) point out that if the metal- 
licity is high, the temperature must also be elevated. 
The ^ 1 0,000 K temperatur e suggested b vlGarcfa Muhoz 
(2007), iMurrav-Clav et al.l (|2009f h and IKoskinen etaL 
(|2010l ) for the thermosphere may not explain the large 
abundances of Si 3+ needed for the observed absorp- 
tion. The ionization energy, AE/k for Si 2+ — >• Si 3+ = 
3.9 x 10 5 K, suggests that the Si 3+ may be produced in a 
s hock between the st ellar and planetary winds. 

iLinskv et all (2010) find no significant Si iv absorption 
by HD209458b. Our best fit to the STIS data predicts 
the average transit depth to be 9^% for the same phases 
as their observations, but ILinskv et al.l (|2010l ), with the 
Cosmic Origins Spectrograph (COS), found 0.2 ± 1.4%. 
These results disagree at the ^1.7 a level and variability 
in the planet's atmosphere may account for the discrep- 
ancy. ILinskv et al.l ((2010) did observe some weak Siiv 
absorption features found at +20 and +40km/s in their 
spectrum, which indicates that some Si 3+ ions remain in 
t he planet's atmosphe re or winds. 

ILinskv et al.l (|2010l ) suggest that the amount of an- 
other Silicon ion, Si 2+ , may vary appreciably over short 
timescales because of changes in stellar wind speed, plan- 
etary mass- loss rate or temperature fluctuations. Con- 
versely to Siiv, for which 2003 STIS observations in- 
dicate strong absorption and 2009 COS observations 
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indicate weak absorption, the Sim absorption is seen 
strongly in absorption in the 2009 COS data and weakly 
in the 2003 STIS data. 

Additional observations will help confirm or refute the 
detection of Si 3+ in the thermosphere of HD209458b. 
The Cosmic Origins Spectrograph (COS) is an ideal in- 
strument with 2 to 10 times the sensitivity of previ- 
ous ultraviolet sp ectrographs (Fronin g fc Greenl [2009). 
COS was used by iLinskv et al.l (j201Q[ ). but only when 
HD209458b was close to a phase of 0.0, 0.25, 0.5 and 
0.75 and once when the planet was at the stellar limb. 
Additional observation at the host's limb would be op- 
timal for Si iv and other optically thin emission lines. 
With enough signal to noise, the light curve may reveal 
asymmetries in the transit having to do with an asym- 
metric spatial distribution of the UV-absorbing cloud. 
The advantage of the limb brightened emission lines is 
that they come from a smaller spatial area of the star 



and therefore probe the spatial distribution of the plan- 
etary atmosphere better than optically thick emission. 
Accurate time resolved transit data may also reveal dif- 
ferences in thermal properties of t he leading and t railin g 
sides of the planet as predicted bv lFortnev et al.l ((2010). 
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